The effects of interactions and disorder in the two-dimensional chiral metal 
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We study the two-dimensional chiral metal, which is formed at the surface of a layered three- 
dimensional system exhibiting the integer quantum Hall effect by hybridization of the edge states 
associated with each layer of the sample. We investigate mesoscopic fluctuations, dynamical screen- 
ing and inelastic scattering in the chiral metal, focussing particularly on fluctuations of conductance, 
Sg(B), with magnetic field, B. The correlation function (Sg(B)Sg(B + SB)) provides information 
on the inelastic scattering rate, , through both the variance of fluctuations and the range of 
correlations in SB. We calculate this correlation function for samples which are not fully phase 
coherent. Two regimes of behaviour exist, according to whether is smaller or larger than tJ_ , 
the rate for inter-edge tunneling, and we give results in both regimes. We also investigate dynamical 
screening of Coulomb interactions in the chiral metal and calculate the contribution to r^ 1 from 
electron-electron scattering, finding oc T 3 / 2 for r^ 1 < t^ 1 at temperature T. 

PACS numbers: 73.40.Hm, 73.23.-b, 73.20.-r, 72.15.Lh 



I. INTRODUCTION 

An unusual type of two-dimensional conductor, the 
chiral metal, may be formed at the surface of a lay- 
ered, three-dimensional system in the presence of a mag- 
netic field which has a component perpendicular to the 
layersErcl. It occurs in a system of weakly-coupled lay- 
ers at magnetic field-strengths for which an isolated layer 
would have an integer quantised Hall conductance. Un- 
der these circumstances, bulk states are localised and the 
extended, chiral edge states associated with each layer 
hybridise to form a two-dimensional conductor at the sur- 
face of the sample. Properties of this chiral metal may be 
probed by conductance measurements with current flow 
normal to the layers. 

Earlyj-experimental work on semiconductor multilayer 
samplesQ established that the integer quantum Hall ef- 
fect of a single layer is indeed robust against weak in- 
terlayer tunneling, and suggested that, within a quan- 
tum Hall plateau, interlayer conductivity in the bulk van- 
ishes in the low»temperature limit. Subsequent theoreti- 
cal discussionsS'B drew attention to the unusual nature of 
surface states in multilayer quantum Hall systems, and 
some aspects of the theory of the chiraLmetal have since 
been explored in considerable detai]d1l3. Recent exper- 
iments by Druist et al have isolated, the contribution of 
surface states to vertical transport^ and demonstrated 
that mesoscopic conductance fluctuations offer. &, way to 
study inelastic scattering in the chiral metaO'tS. Con- 
duction by surface states has also been invoked in the 
interpretation of further experiments on semiconductor 
multilayer samples^, and on tie bulk quantum Hall ef- 
fect in both organic conduciarsEj and an inorganic quasi- 
two dimensional conductorO. 

In this paper we extend the theoretical treatment of 
mesoscopic chiral conductors in two respects. First, we 
develop some aspects of the theory of conductance fluctu- 
ations which we hope will be useful for the interpretation 



of experiments. Second, we study electron-electron inter- 
actions in the dirty chiral metal. 

The existing understanding of conductance flue 
tions in the chiral metal is based on calculations! 
for the variance of sample-to-sample conductance fluctu- 
ations in systems that are fully phase-coherent. From 
the results of those calculations, estimates have been 
madeO for the variance of conductance fluctuations mea- 
sured as a function of magnetic field in samples that are 
only partially phase-coherent, and also for the range in 
magnetic field of conductance correlations. These es- 
timates determine the relevant scales but do not yield 
numerical coefficients or functional forms. Moreover, 
they depend on a continuum treatment of the system 
of coupled edge states, and apply if the inelastic scat- 
tering rate, is smaller than the inter-edge tunneling 
rate, t± , while existing experiments appear to be in 
the opposite regimeEZHl3. We supplement this past work 
by calculating in full the conductance correlation func- 
tion (5g(B)5g(B + SB)), both for r^ 1 < tJ 1 and for 

Separately, we consider conductance as a function of 
Fermi energy. For conventional disordered conductors, 
changes in Fermi energy provide a way to sample the 
conductance distribution. In contrast, we find that the 
conductance of a partially-coherent chiral metal sam- 
ple, though dependent on the disorder realisation, is un- 
changed by small variations in Fermi energy. One conse- 
quence is that thermal smearing of the electron distribu- 
tion function does not decrease the amplitude of conduc- 
tance fluctuations in the chiral metal. 

Turning to electron-electron interactions, we estab- 
lish the form of the dynamically screened interaction 
and examine some of its consequences. Most impor- 
tantly, we calculate the inelastic scattering rate due to 
electron-electron interactions, obtaining cx T 3 / 2 for 
Tin 3> t±. Wc also examine other interaction effects 
which are known to be of interest in conventional disor- 
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dered conductors^. Interactions in diffusive metals are 
responsible for a zero-bias anomaly in the tunneling den- 
sity of states^, but - as noted previouslya, and as we con- 
firm - this is absent from the chiral metal. Similarly, the 
leading interaction contribution to the mean conductance 
has a singular temperature dependence in conventional 
disordered conductorsE-aES, but we find that it cancels in 
the chiral metal. 

The remaining sections of this paper are organised as 
follows. We review briefly past work and collect our most 
important new results in Sec. II. We describe calcula- 
tions on conductance fluctuations in Sec. Ill, and treat 
electron-electron interactions in Sec. IV. We summarise 
and discuss recent experiments in Sec. V. 



II. BACKGROUND AND MAIN RESULTS 



It is useful first to introduce some notation. Consider 
a layered sample which exhibits the bulk quantum Hall 
effect, with one Landau level below the Fermi energy. A 
chiral metal is formed at the sample surface, from hy- 
bridisation of one edge state per layer. We denote the 
interlayer spacing by a, the chiral velocity by v and the 
interlayer diffusion constant by D. The electron density 
of states per unit area is n — l/hva and the surface 
conductivity in the transverse direction, as given by the 
Einstein relation, is a = D/va, in units of e 2 /h, which 
we use throughout this paper. 

Phase coherent mesoscopic conductors are known to 
have three distinct mesoscopic regimes, according to rela- 
tive sample dimensional^. For definiteness, suppose that 
the sample is a cylinder of height L and circumference 
C, with its axis perpendicular to the layers. Take the 
number of layers to be N, so that L — Na. Two in- 
trinsic length scales arise, with which L should be com- 
pared. The first of these involves the mean conductance 
of the sample, which if Ohm's law applies is (g) = Ca/L: 
the quasi one-dimensional localisation length, £ = 2Cer, 
is the value of L at which (g) is of order the quan- 
tum unit of conductance. Samples with L 3> £ have 
states localised in the vertical direction and are quasi one- 
dimensional insulators. Samples with L <C £ are metallic 
and have Ohmic dependence of (g) on sample dimensions. 
A second length scale influences the amplitude of con- 
ductance fluctuations. It is the distance L\ that an elec- 
tron diffuses in the transverse direction during the time 
taken to circumnavigate the sample: L\ — (DC/v) 1 / 2 = 
(acrC) 1 ' 2 . Samples with L\ <C L -C £ are quasi one- 
dimensional metals without time-reversal symmetry and, 
in common with other examples of such systems, have a 
variance for the conductance of ((Sg) 2 ) = 1/15. Sam- 
ples with L < Li are two-dimensional chiral metals. If 
contacts are attached to these samples for conductance 
measurements, electrons will escape by diffusion to the 
contacts before completing a circuit of the sample in the 
chiral direction. As a result, multiple scattering of an 



electron from a given impurity is completely supressed. 
Such samples present a number of essentially indepen- 
dent parallel strips for conduction in the transverse di- 
rection, each having a width in the chiral direction given 
by the distance an electron propagates before reaching a 
contact. This width is L 2 v/D = C{L/Li) 2 . The conduc- 
tance of the two-dimensional chiral metal may be thought 
oiQ as a sum of contribiitioiis from each such strip, and 
its variance has the value&E£l (Sg 2 ) = L\jJ>L 2 . 

Our concern in the following is with samples that are 
not completely phase coherent. They have two differ- 
ent metallic regimes, analogous to those described above 
for phase coherent conductors. We restrict our attention 
to one of these, the incoherent two-dimensional metal, in 
which T in < C/v so that dephasing occurs before circum- 
navigation: this is the simpler regime theoretically and 
seems the one relevant for existing experiments. Cho, 
Balents and Fishertil have estimated the size of conduc- 
tance fluctuations in the incoherent metal by viewing 
each phase-coherent region as a classical resistor and the 
whole sample as a random resistor network. We sum- 
marise their argument here for convenience. Let the 
number of regions into which the sample is divided be 
n x in the chiral direction, and n z in the transverse direc- 
tion. Let the conductance of a single region have mean 
value ffpatch, with fluctuations of magnitude <5g pa tch; pro- 
vided (Jgpatch *C gpatch, the corresponding resistances 

are -R pa tch = l/ffpatch and <> -Rpatch = Rpatch/ <7 pa t c h- 

A strip of n z such regions forms a single conduction 
path in the transverse direction, having mean resistance 
-Rstrip = n z R patch with fluctuations SR 2 trip = n z SR 2 &tch ; 
the corresponding conductances are g s t r i P = 1/-Rstrip 
and <5<? s trip = SR stl i P / R 2 tl[p . Since the sample consists 
of n x such paths in parallel, its conductance has mean 
value (g) = n x g strip = (n x /n z )g pgbtch and fluctuations 
(Sg 2 ) = n x Sg 2 trip = {n x /nl)8gl &tch . The implications 
of this final result depend on the size of a patch in the 
transverse direction, compared to the layer spacing, a. 

If <C tJ 1 so that the phase coherence length in the 
transverse direction is much larger than the—layer spac- 
ing, which is the case considered previously^, then one 
expects a continuum treatment of the system to be ad- 
equate. In this case, n 
in addition, n x = C/vt, 
and introducing the inelastic scattering length in the chi- 



- - L 2 /Dr in and 5g 2 p&tch ~ 1; 
Combining these expressions 



ral direction, L 



one arrives at the conclusion 



(Eq. (4.1) of Ref. [Tl|) that 



(Sg 2 



c 



1/2 



N 



3/2 



(2.1) 



Experimental studies use a magnetic field with a compo- 
nent normal to the sample surface to.|generate conduc- 
tance fluctuations in the chiral metal3'EZI; the field scale 
of conductance correlations is set by the flux quantum, 
$0, and the area of a coherent region, being (Eq. (4.4) of 
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(2.2) 



The resu lts we obtain in Sec. Ill substitute for Eqns. (2.1) 
and (2.2) the correlation function 



(2.3) 



where the function f(y) is 

/(») = ir- 1 ' 2 I + (2.-II 

with the scaling variable 



7T 2 ;3 /&B 

12 m V $o 



(2.5) 



The opposite limit of weakly coupled edges, r in ^> 
tJJ , has not previously been examined theoretically but 
seems to be important experimentallyllSEa. Adapting the 
argument outlined above, we take, as elements of a clas- 
sical resistor network, regions again of length l in in the 
chiral direction but now of width a in the transverse di- 
rection. Then n x — C/li n and n z — N. The mean 
conductance of one such patch is small in this limit. Wc 
expect fluctuations to be of the same order, and take 
^Spatch ~ ffpatch = ((g)n z /n x ) 2 . In this way we obtain in 
place of Eq. (O) 



(Sg 2 



N ' C 



(2.6) 



The field scale of conductance correlations is again set by 
the flux quantum and the area of a coherent region, but 
Eq. 2.2 is replaced in this limit by 



SB ~ $ /('m«) • 



(2.7) 



We calculate in Sec. Ill the conductance correlation func- 
tion for weakly coupled edges, finding in agreement with 
these estimates 



(Sg(B ± )Sg(B ± + SB)) 



NC 1 



(2- 



where z = 2irSBli n a/^Q. 

The inelastic scattering length, li n , appears as a phe- 
nomenological parameter in the expressions given above 
for the correlation function of conductance fluctuations. 
The contribution to inelastic scattering from electron- 
electron interactions is known in non-chiral dirty metals 
to have a characteristic temperature dependence, reflect- 
ing the form of the dynamioaJly-,S|Gseened Coulomb inter- 
action in a diffusive systen£3'EZrE3, and it is interesting 



to study inelastic scattering microscopically for the chi- 
ral metal. We do this in Sec. IV, examining dynamical 
screening and using the results to calculate the inelastic 
scattering rate, r^I . We obtain in the regime r^ 1 <C rj 1 



' D 1 / 2 



fc R T 



3/2 



(2.9) 



where c w 1.5 is a dimensionless coefficient. A sim- 
ple interpretation of this result can be given, following 
similas-arguments established in the theory of diffusive 
metalsc3. As a starting point, note that (in both chiral 
and non-chiral diffusive metals) the dynamically screened 
interaction strength at long wavelengths is independent 
of the bare interaction strength. In consequence, the 
inelastic scattering rate depends only on properties of 
the non-interacting system and temperature: fvr~^ is of 
order the single-particle energy level spacing in a sys- 
tem with dimensions determined by the distance an elec- 
tron moves in the time h/k^T. These dimensions are 
L X (T) ~ hv/k B T and L Z (T) ~ (hD/hsT) 1 / 2 in the chiral 
and transverse directions respectively, and the estimate 
t^ 1 - {hnL^L^T))- 1 ) is consistent with Eq. (§J). 
An equivalent argument applied to non-chiral, diffusive 
conductors in d-dimensions yields t^ 1 - T d / 2 , which 
is corrected in two dimensions to r-" 1 ~ T|log(T)| by 
detailed calculationsE3Ea. The factor of log(T) in this 
last expression reflects divergent contributions to scat- 
tering in the diffusive two-dimensional metal, from pro- 
cesses with small energy transfer. The same scatter- 
ing processes are responsible for a difference, in diffusive 
two-dimensional metals, between the temperature depen- 
dence of the inelastic (or out-scattering) rate, which acts 
as a cut-off for conductance fluctuations, and the dephas- 
ing rate, which acts as a cut-off for the weak localisation 
effects that result in negative magnetoresistance. In the 
chiral metal, as for conventional, diffusive metals in more 
than two dimensions, small energy transfer processes are 
not dominant and there is a single relaxation rate. 

Inelastic scattering in weakly coupled edge states has 
been discussed in Ref. ||, where it is noted that one ex- 
pects Tj" 1 °t T from perturbation theory in the interac- 
tion stength, applied to a single edge. 



III. CONDUCTANCE FLUCTUATIONS 
A. Model 

We study conductance fluctuations in the chiral metal 
using a single-particle description of the system, supple- 
mented by a relaxation rate to represent inelastic scat- 
tering. We treat a sample with layers lying in the x — y 
plane and consider a surface in the x — z plane. Electrons 
on the surface then have a continuous coordinate, x, in 
the chiral direction parallel to the edges of layers, and 
an integer coordinate, n, labeling layers in the transverse 
(z) direction, which we combine as r = (x,n). We take 
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the interlayer tunneling energy to be t and represent a 
magnetic field component, B±, normal to the surface us- 
ing the vector potential A = B^ank. The Hamiltonian 
H acts on a wavefunction tp n (x) according toBtJ 

(Hip)n(x) = v(-ihd x + eB±an)ip n (x) 
-t[lp n +l(x) + 1p n -t(xj\ + V n (x)lp n (x) , (3.1) 

(3.2) 

where V n (x) is a random potential arising from impurities 
and surface roughness. We choose V n (x) to be Gaus- 
sian distributed with short-range correlations, so that 
(V n (x)) = and {V n (x)V m (x')) = AS nm 6(x - x 1 ). 

Denoting the Green's function for Eq. ( |3.l| ) by g = 
(z — -ff) -1 , we require its disorder average 

G(z;n,r 2 ) = (g(z;n,r 2 )) , (3.3) 

the diffusion propagator 

K{w\ ri, r 2 ) ee (g{u> + iO; n, r 2 )ff(-«0; r 2 , r x )) , 



(3.4) 



and its Fourier tranform 



dx 



(3.5) 

These quantities are known from calculationsBB using 
the usual expansion in powers of the Green's function 
for the disorder-free system. Without disorder or a 
normal magnetic field component, eigenstates are plane 
waves and the electron dispersion relation is e(k) = 
hvk x — 2tcos{k z a). This dispersion along with the chi- 
ral motion leads to an open Fermi surface on which all 
electrons have the same x-component of velocity. In the 
presence of a normal magnetic field component of dimen- 
sionless strength b = eB±va/t, the eigenfunctions ipn{x) 
in the pure system are 



ip n (x) 



1 



cj> a (n) , 



(3.6) 



where a is an integer and 4> a (n) — J„_ Q (2/6), the Bessel 
function of order n — a. The associated energy eigenval- 
ues are hvk + e a with e a = abt. 

The treatment of disorder is simple in some respects 
because chiral motion prevents multiple scattering of an 
electron from any particular impurity. As a result, the 
average one-particle Greeks function is given exactly by, 
the Born approximatiorJ3ll3. In this way, one arrives att3 



G R (E;r 1 ,r 2 ) = 

2tt J-oo ^ E + iA/(2Hv) - (hvk + e a ) 



(3.7) 



where we use the standard notation, G R (E) = G(E + iO) 
and G A (E) = G(E - iO). From the value of the self- 
energy, the elastic scattering time is r e ; = h 2 v/A, and the 



elastic scattering length is l e i = (hv) 2 /A. In the absence 
of a normal magnetic field component, G R (E; r 1? r 2 ) is 
translationally invariant and has the Fourier transform 



G R (E; k) = [E + t'A/(2fi«) - (hvk x 



2tcos(k z a)] 1 . 

(3.8) 



For non-zero B±, G R (E; ri, r 2 ) acquires a phase under 
translations, which is central to calculations of conduc- 
tance fluctuations with magnetic field: 



G R (E; 0, r) = exp(i(3an'x)G R (E; r',r + r'). 



(3.9) 



where (3 = 2irB±/&o, r = (x,n) and r' = (x',n'). 

The diffusion propagator is similarly given by a a sum 
of ladder diagrams, and has a form that reflects the com- 
bination of chiral motion in x and diffusion in 2. At small 
wavevectors it is 



K(u,k) = {hvy^hDk 2 - i{uj + hvkx)]' 



(3.10) 



The diffusion constant has the value D — 2(at) 2 v/A 
in the absence of a normal magnetic field component, 
and has the field dependence^ D(B±) = D(0)/[1 + 
(B±/ Bq) 2 ], where the field scale for magnetoresistance 
is B = $ /al e i. 

Because of the simplifications in the treatment of dis- 
order that follow from chiral motion, the results above 
hold both at weak disorder (t » A/(hv)) and at strong 
disorder (t <C A/(hv)). In a semiclassical description, 
electron motion follows a random walk in the transverse 
direction. At weak disorder, steps of the walk have speed 
at/H and duration r e ;; the rate for inter-edge tunneling 
is rj 1 = t/h. At strong disorder, steps are of length a 
and duration t±; from the value of D, we deduce that 
the rate for inter-edge tunneling is then tJ 1 = (t/h) 2 T e [. 



B. Conductivity 

The component of the current density operator in the 
transverse direction is 



Cdt 

jz{r) = —[i)l(x)ip n+1 (x) - ij}l +1 (x)ip n (x)] 



(3.11) 



The Kubo formula for the transverse conductivity in a 
system of area A is 



(3.12) 



Evaluating the disorder average of this expression, one 
finds 



e 2 D 
h va 

as expected from the Einstein relation. 



(3.13) 
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C. Conductance fluctuations in magnetic field 



1. Strongly coupled edges: 



ConducteuiGe fluctuations have been studied 
previouslyfl'Lni3 for the two-dimensional chiral metal in 
the absence of inelastic scattering, in most detail by 
Gruzberg, Read and Sachdevt3, who discussed the vari- 
ance and its dependence on sample geometry. Here we 
present calculations which include inelastic scattering, 
obtaining the correlation function 



F(SB) = (Sg(B)Sg(B + SB)) 



(3.14) 




4 3 

FIG. 1. The Feynman diagram which is relevant for the 
conductance fluctuations: current operators are represented 
by circles, one-particle Green's functions by straight lines, 
and diffusion propagators by wavy lines. 

As a result of chiral motion, and provided that elec- 
trons do not circumnavigate the sample phase coherently 
{C/v > Ti n or C/v > L 2 /D), the only one-loop diagram 
contributing to F{5B) is the one shown in Fig. [|. It is 
most conveniently evaluated in real space, and represents 
the expression 

F(SB) = ( Y,K S B{ri,V2)K* SB (r A ,rz) 

rii 

J{r l7 r 4 )J{r 2 ,r 3 )d{xi} . (3.15) 

Here, Kg B (ri, ^2) stands for the zero- frequency diffusion 
propagator, generalised to the situation of interest, in 
which the two Green's functions entering it are evaluated 
for different normal magnetic field strengths, 



K SB (r,r') = {gX(r,r')gA(r',r)), 



(3.16) 



and in which inelastic scattering is included. The other 
factors arise from a combination of single-particle Green's 
functions and current operators, and are short range: 
J(r,r') = (hA 2 /2TTL 2 )\C(r,r')\ 2 , with 



C(r,r') = J J2G A (r,r 1 ) Jz (r 1 )G R (r i y)dx 1 



(3.17) 



Our approach to evaluating F(5B) is different in the two 
regimes, <C 7"j_ and t±, which we consider 

separately. 



In this regime, the discreteness of the system in the 
transverse direction may be neglected. It i s suff icient to 
approximate the short-range terms in Eq. ( 3.15 ) by 



J(r,r') = J 5(r-r' 



where 



Jo = 



e 2 Ajatf 
h AL 2 



(hv) 2 



(3.18) 



(3.19) 



The diffusion propagator, generalised to include the mag- 
netic field difference SB and inelastic scattering, is most 
easily calculated within a continuum treatment of the 
transverse direction by solving, for x > 0, the differential 
equation 

vd x K SB {0, r) = [Dd 2 + ivSPz - r^KssiO, r) 

(3.20) 

with the initial condition Kg B {Q]x, z) = (hv)~ 2 5{z) at 
x = 0. Here, the inelastic scattering rate, r" 1 , has been 
introduced, 

and the magnetic field difference enters through 8j3 = 
(2it5B/&o). To describe a sample connected to contacts, 
absorbing boundary conditions (d z KsB{0, r) = 0) should 
be imposed at z — and z = L. However, provided the 
size of a phase-coherent region is much smaller than the 
sample size, so tha t DT j n <C L 2 , it is sufficient to use 
the solution of Eq. ( 3.2d| ) for a system of infinite extent 
in t he z-d irection when evaluating the right-hand side of 
Eq. (|3.15|). This solution is 



K SB (0,r) = {hv) 



AttDx 



1/2 



exp(-S) 



(3.21) 



where 



S 



4Dx 



x 

hn 



iSPxz D(SP) 2 x 3 



12v 



(3.22) 



Combining Eqns. ( |3.15| ), ( |3.18| ), and ( |3.2l| ), we obtain 

/>oo />oo 

F(6B) = J$LC dx dz\K SB (0,r)\ 2 

Jo J-oc (3.23) 

and so find the correlation function for conductance fluc- 



tuations given in Eq.(2.3) 

The form of the scaling function f(y) is shown in 
Fig. H It is normalised so that /(0) = 1, and decreases as 
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/ oc y for y >> 1, so that the conductance correla- 
tion function decays as F(SB) oc j^i?!" 1 / 3 for large \SB\. 



D. Conductance fluctuations in energy 
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FIG. 2. The scaling function f(y) 



It is also of interest to consider the dependence of con- 
ductance fluctuations on Fermi energy. In particular, the 
range of conductance correlations in energy for a system 
at zero temperature will determine the extent to which 
the amplitude of fluctuations measured at finite temper- 
ature is reduced by thermal smearing of the electron dis- 
tribution. We find that conductance fluctuations are per- 
fectly correlated in energy. To show this, we return to 
Eq. (3.15) but replace the magnetic field difference SB 
with an energy difference SE. The diffusion propagator 
at finite energy difference simply acquires a phase: 



K se (t,t') = K(r,r')e 



i5E(x'-x)/(hv) 



(3.28) 



2. Weakly coupled edges: 



>t7 



If adjacent edges are sufficiently weakly coupled by 
tunneling, there is mesoscopic regime in which the inelas- 
tic scattering rate is small, in the sense that r^ 1 <C tT, , 
but the interlayer tunneling rate is even smaller. It is 
sufficient in this regime to calculate the conductance cor- 
relation function at leading order in the interlayer cou- 
pling, t. Doing so, we retain t in the current operators 
but use Green's functions for a system of isolated edges. 
The disorder-averaged single-particle Green's function in 
this limit 



is 



G R (E-v 1 ,v 2 ) 
(ihv)- 1 exp(-[(2Z 



0(x 2 - xi)S nitn2 x 
_1 - i(3ani}[x2 - x\]) , 



(3.24) 



where Q(x) is the step function. The short-range terms 
associated with current operators in Eq. (3.15) may then 
be approximated as 



J{v,v') = ^5{x 



X')[S.. 



n,n' + l 



(3.25) 



The diffusion propagator reduces for uncoupled edges to a 
function which has a phase determined by SB and which 
decays in amplitude with distance only as a result of in- 
elastic scattering. We take it to be 

K SB (ri, r 2 ) = 9(x 2 - xi)5 nuri2 x 

{hv)~ 2 expf-pr 1 - i8Pam][x 3 - a*]) . (3.26) 

Combining Eqns. ( pTl5j ), ( |3~25| ), and ( gjgg ), we obtain 



F(SB) = -J 2 NC x 
dx[K S B(0,n;x,n)Kg B (0,n + l;x,n + 1) + c.c] (3.27) 



and hence arrive at the correla tion function for conduc- 
tance fluctuations given in Eq.(|2.8|). 



This phase cancels when Kse( t i> r 2 ) and K$ E (r4, r 3 ) are 
combined, since the factors J(r,r') associated with the 
current operat ors ar e short-range, so that ri « r4 and 
r 2 w r 3 in Eq. ( 3.15 ). As a consequence, the correlation 
of conductance fluctuations is independent of SE. The 
physical reason for th is b ehaviour is that, if ij) n {x) is an 
eigenfunction of Eq. (3.1) for some realisation of disor- 
der, with energy E, then ip' n (x) = e lux ' v ip n (x) is also an 
eigenfunction, with energy E' = E + Huj, provided uoC/v 
is an integer multiple of 27r, so that periodic boundary 
conditions in the chiral direction are satisfied. Hence, 
apart from the phase factor, which does not affect the 
conductance, states are perfectly correlated in energy. As 
a result, the only consequence of a change in temperature 
is a change in the inelastic scattering rate. 



IV. EFFECTS OF INTERACTIONS 



A. Polarization and screening 



In this section we study several aspects of electron- 
electron interactions in the chiral metal. As a first step, 
it is necessary to examine screening of the Coulomb in- 
teraction between a pair of electrons, by other electrons 
in the Fermi sea. One can expect the frequency and 
wavevector dependence of the screened interaction to be 
significant. In conventional conductors, the consequence 
of disorder and the resulting diffusive motion of charge 
(with diffusion constant D) is that screening of a poten- 
tial having wavevector q is suppressed^ at frequencies 
higher than Dq 2 . And in a one-dimensional chiral metal 
at the edge of a two-dimensional integer quantum Hall 
system, electron-electron interactions modify the disper- 
sion relation for edge magneto-plasmonstJ. We find in 



G 



the two-dimensional chiral metal that dynamical screen- 
ing at finite wavevector interpolates, as a function of the 
direction of the wavector, between these two types of be- 
haviour. 

Our results may be obtained either by using the ran- 
dom phase approximation within a diagrammatic calcu- 
lation, or more transparently from a hydrodynamic ap- 
proach, as follows. 

Consider the response of the chiral metal to a space- 
and time-dependent external potential. Let p(r, t) be the 
screening charge density induced in the presence of an 
electric field that has components (E X ,E Z ) within the 
surface. We require an equation of motion for p(r,t), 
which we derive in the usual way, by considering J(r,i), 
the deviation from equilibrium of the current density, and 
using the continuity equation. Treating the transverse 
direction as continuous, one has 



and 



J x {r,t) = vp(r,t) 



J z (r,t) = -Dd z p(r,t) + aE z 



(4.1) 



(4.2) 



so that the effect of the transverse component of the elec- 
tric field is simply to generate a transverse current den- 
sity. The effect of the component of the electric field in 
the chiral direction is rather different: referring to the 
three-dimensional sample as a whole, E x generates Hall 
currents within each layer, which transport charge be- 
tween the bulk and the surface. The resulting change in 
surface charge density may be represented by adding a 
term to the continuity equation, so that it becomes 

d t p(r,t) = -V-3(r > t) + —E x , (4.3) 

a 

where uh is the quantised Hall conductance of a layer, 
and (as in preceding sections) a is the layer sp acin g. It 
i s he lpful to express the conductivities in Eqns. (4.2) and 



( |4.3|) as a = e nD and an = e /h, and to recall that 
the de nsit y of st ates i s n = (hva)^ 1 . Then, combining 
Eqns. (4T), (^2) and (|f ,3|), the charge density evolves as 



dtp(r, t) = 

-vd x p{r 7 t)+Dd 2 zP (r,t) 



2 n[vE x — Dd z E z ]. (4.4) 



The electric field appearing in these equations arises 
from a scalar potential $ tot (r,i) which is the sum of an 
externally imposed potential, $ ext (r, t), and the poten- 
tial <& sc ,.(r,i) due to the screening charge, p(r,t). These 
last two quantities are related in the standard fashionE3 
by the Poisson equation for a three-dimensional electro- 
static problem, with coordinates r = (x, z) and y, in 
which the charge density is p(r,t)S(y) and <& scr (r,i) is 
found by evaluating the three-dimensional potential at 
y = 0. Taking Fourier transforms, defined according to 



p(u,k) = J d 2 r J dte^ k r+wt V(r, t) , 



(4.5) 



and similarly for other quantities, one finds 

<S> scr (iu, k) = l7o(k)p(w, k)/e 2 , (4.6) 



where 



Co(k) 



2e r e |k| 



(4.7) 



is the Fourier transform of the unscreened interaction po- 
tential betwen a pair of electrons. Expressing the electric 
field in term s of $ tot (r,£) and solving the Fourier trans- 



form of Eq. (4.4), we find 



*w(u;,k) 



$ezt(w,k) 



i + [/ (k)n( w ,k) 



where 



, \ Dk z - ivk x 
i-L[oj, k) = n o 

Dk z — iui — ivk x 



(4.8) 



(4.9) 



This result simplifies in several limits. Static screening 
(u> = 0) is isotropic and as given by Thomas-Fermi the- 
ory, with an inverse screening length n = e 2 n/2e r eo. Dy- 
namical screening of a potential with Fourier components 
only in the transverse direction (k x — 0) is exactly as in a 
non-chiral disordered, two-dimensional metal. Response 
to a potential with Fourier components only in the chiral 
direction (k z = 0) is undamped, and excitations have the 
dispersion relation: lu = —v(k x + k). 

In the following, we shall treat interactions in the chi- 
ral metal using the Matsubara formalism. Within this 
approach, we find for the polarisation operator 



n(ttj n ,k) = n 



ivk x sgn(u n ) 



Dk z 



+ |w„| - ivk x sgn(uj n ) (4.10) 



with u) n = 27mksT, which has, as its analytical contin- 
uation, Eq. ( [4.91 ). Similarly, the screened Coulomb inter- 
action potential is 

k) = [/o(k)/[i + [/ (k)n^„,k)]. 



B. Tunneling density of states 

The enhanced interaction between slowly relaxing den- 
sity fluctuations in a conventional, diffusive metal is re- 
sponsible for_a. zero-bias anomaly in_the tunneling den- 
sity of statea^j. Balents and Fishein have argued that 
such a zero-bias anomaly should not be expected in the 
chiral metal. Physically, charge that tunnels into the 
system is swept away at the chiral velocity; mathemat- 
ically, the terms in perturbation theory associated with 
this anomaly yield only smooth contributions for the chi- 
ral metal. 
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In the present subsection we examine in more detail the 
reasons for this. It is sufficient for the purpose to consider 
a static interaction, U (r), with finite range, between elec- 
trons separated by a distance r, and to calculate at ficsjt 
order in U(r) the disorder-averaged exchange energyEl 
Let E(-E) be the disorder-average of this quantity for-aii 
electron at energy E interacting with a filled Fermi seaE3. 
The change in the density of states due to interactions is 



6n(E) 



dE(E) 
1 dE 



(4.11) 



Viewed in this way, the zero-bias anomaly in conven- 
tional, diffusive metals arises because states lying close 
in energy have wavefunctions which are highly correlated 
in space, enhancing interactions. Moverover, the wave- 
function correlations depend strongly on the energy sep- 
aration of the states, so that T,(E) varies rapidly with E 
if E is close to the Fermi energy, Ep. 

Let { x f' Q ,(r)} and {E a } be single-particle eigenfunc- 
tions and eigenenergies in a disordered chiral metal. The 
disorder-averaged exchange energy is 

-i />Ep — E /> 

£(£) = / du / d 2 rU(r)A(u, r) , 

nj-oz J (4.12) 

where the two-particle spectral function, A(cj,t), is de- 
fined by 

A(u, r) = (J2 S(E q )6(cj - ^)*^(0)* a (r)^(r)^(0)) . 

a/3 

This spectral function can be calcul ated from the real 
part of the diffusion propagator, Eq. (3.4), or by consid- 
ering the time-Fourier tranform of the density associated 
with a spreading wavepacket. It has the form 



1 / 



1/2 



exp(- 



4Dx 



'-) 



(4.13) 



a(hv) 2 V47rD|a;| 

Combining Eqns. pll| ), ( glf ) and Cl| ), one sees that, 
if the potential has range R, then Sn(E) is independent 
of E for \E- E F \< hv/R. 

Inspecting Eq. ( 4.13| ), it is clear that eigenfunctions 
in the chiral metal do in fact have strong correlations: 
A(w,r) with r = (x,z) diverges as \x\ — > for z = 0. 
These correlations, however, depend only weakly on en- 
ergy, and so do not produce a large change in the tun- 
neling density of states. 



C. Inelastic scattering rate 



The following subsection is devoted to a calculation 
of the contribution from electron-electron interactions to 
the inelastic scattering rate, r~ 1 , which acts as a cut-off 
for conductance fluctuations. 

We compute r^ 1 by finding at leading order the in- 
teraction contribution to the irreducible vertex that ap- 
pears in the ladder sum for the diffusion pcanaffator used 
in calculations of conductance fluctuationsEZrE^I. The two 
single-particle (advanced and retarded) Green's functions 
from which this diffusion propagator is constructed are 
associated with distinct measurements of the conduc- 
tance. In consequence, one should not include in the 
irreducible vertex diagrams in which the interaction line 
connects the two single-particle propagators. The rel- 
evant diagramsE3 are shown in Fig. ^, in which the full 
lines represent single-particle propagators, wavy lines are 
the screened Coulomb interaction, single dashed lines are 
impurity averages, and double dashed lines are diffusion 
ladders. We evaluate these diagrams treating the trans- 
verse direction as continuous, so that our results are re- 
stricted to the regime in which edges are strongly cou- 
pled: rj 1 ^s> Tj„ . We set the energy and wavevector 
differences, lj and q, at the vertex to zero, express sums 
on Matsubara frequencies as contour integrals, and com- 
bine terms. We expect that the dominant contribution to 
the resulting integral will be from small wavevectors and 
energies, and use the asymptotic form of the integrand 
appropriate for this regime, obtaining 



d 2 q 



dx x 



irhn J (2tt) 2 
[coth(a;/2fcBT) - tanh((a;)/2fcBT)] x 

x(hDq z 2 ) 2 

[(as + Hvq x ) 2 + (hDq z 2 ) 2 ][(hvq x ) 2 + {KDq 2 ) 2 ] ' 



(4.14) 



The dependence of this expression on T, D and v is 
made apparent by introducing the scaled variables: X = 
x/k B T, Q x = hvq x /k B T and Q z = (HD /k B T) l / 2 q z . The 
integral can be evaluated n ume rically, yielding the ex- 
pression for Ti n given in Eq. (2.9). The physical interpre- 
tation of this result has been discussed in Sec. II. 
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e + co — n, 

(i) (ii) 



(iii) 



(iv) 



FIG. 3. The diagrams from which the inelastic scattering rate is calculated. 



D. Interaction corrections to conductivity 

Electron-electron interactions are known in diffusive 
conductors to give rise to temperature-dependent con- 
tributions to the conductivity that are singular in the 
low-temperature limiiOcj The origin of these interac- 
tion corrections to the conductivity is similar to that of 
the zero-bias anomaly in the tunneling density of states. 
Since the chiral metal has no zero-bias anomaly, one ex- 
pects it to have no singular interaction contribution to 
the conductivity. In the following subsection we outline 
calculations that confirm this at leading order. The pre- 
sentation follows-elosely that of Altshuler, Khmel'nitskii, 
Larkin, and Leec3 for the diffusive metal. 

To compute the interaction correction to the conduc- 
tivity of the chiral metal using the Matsubara technique, 
it is convenient to use the Kubo formula in the form 

a = lim — (4.15) 

where A(lu) is the analytic continuation in lu of 



K(iu) n ) = 



The diagrams that contribute to A(iu> n ) at leading or- 
der in the screened interaction strength can be classified 
according to the number of diffusion propagators each 
one contains. Those with the most diffusion propagators 
are most singular in the low-temperature limit. They 
can be separated into two categories. Diagrams in the 
first category are illustrated in Fig. ^, corresponding to 
Figs. 5(a), 5(b) and 5(c) of Ref. |26[ It has been shown 
for a diffusive conductor in Ref. |26| that these cancel, 
and we find the same to be true in the chiral metal. Di- 
agrams in the second category are illustrated in Fig.|^, 
corresponding to Figs. 5(d) and 5(e) of Ref. |6[ For the 
two-dimensional diffusive metal, these diagrams make a 
contribution to the conductivity which is logarithmic in 
temperature. For the chiral metal, their contribution is 
formally proportional to T 1 / 2 , but with a numerical co- 
efficient which we find to be zero after integrating over 
the component in the chiral direction of the momentum 
transferred by the interaction. 



(Tr[j z g(i[e m + u n ])j z g(ie m )]) 



(4.16) 



FIG. 4. Diagrams for the interaction correction to the conductivity which cancel. 




FIG. 5. The remaining leading-order diagrams for the interaction correction to the conductivity. 
V. DISCUSSION ACKNOWLEDGMENTS 



Our results, as set out in Sec. II, are a contribution 
to the understanding of transport in the chiral metal. 
In many senses, it is a particularly simple example of a 
disordered conductor. Quantum interference effects are 
suppressed by chiral motion, and interaction effects are 
less singular than in diffusive two-dimensional conduc- 
tors. We have provided a rather detailed treatment of 
conductance fluctations, which we hope will be tested 
experimentally. We have also identified a distinct meso- 
scopic regime, for weakly coupled edges, in which the 
intcr-cdgc tunneling rate, tJ 1 , is smaller than the inelas- 
tic scattering rate, . In this regime, transport is con- 
trolled by the properties of an isolated, one-dimensional 
chiral metal. 



Recent experimentsli3£3 appear, in fact, to probe be- 
haviour of weakly coupled edges, though their analysis 
has necessarily so far been based on theory for the oppo- 
site regime. According to that analysis, the inelastic scat- 
tering length in the transverse direction, estimated from 
the variance of conductance fluctuations using Eq. (2.1), 
is shorter than the inter-layer spacing, even at the low- 
est temperatures investigated. Moreover, this estimate 
of the inelastic scattering length is inconsistent with the 
value deduced from the cor relation field for condcutancc 
fluctuations, using Eq. ( |2.2[ ), as one might expect if the 
theory used is, in fact, inappropriate. It will be interest- 
ing to find whether these discrepancies are resolved by 
using the theory we have presented to analyse the ex- 
periments. If so, such measurements provide a way to 
investigate interaction effects in a one-dimensional sys- 
tem of chiral fcrmions. 
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